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SUMMARY 


A  research  proposal  to  investigate  the  “Compressibility  Effects  on  Dynamic  Stall 
of  Oscillating  Airfoils”  was  submitted  to  ARO  and  the  project  was  funded  in  April 
1986.  The  aim  was  to  obtain  a  basic  understanding  of  the  effect  of  compressibility 
on  the  phenomenon  of  dynamic  stall  under  typical  flight  conditions  encountered  by  a 
helicopter  in  forward  flight,  so  that  eventually  a  means  for  its  control  can  be  devised 
and  thus,  its  flight  envelope  can  be  expanded. 

The  initial  phase  of  the  study  was  devoted  to  building  a  drive  system  to  produce 
the  necessary  unsteady  airfoil  motion.  A  novel  design  was  arrived  at  after  reviewing 
the  various  possibilities  and  was  built.  It  uses  a  four-bar  chain  mechanis.m  of  which 
the  airfoil  is  one  of  the  links.  The  drive  can  produce  a  sinusoidal  variation  of  the 
angle  of  attack  as 

a  =  0(1  +  amsinujt 

with  the  mean  angle  of  attack  Oco  continuously  variable  from  0  to  16  degrees,  the 
amplitude  of  oscillation  Qm-  from  2"  -  10  ’  and  the  frequency  from  0  -  100  Hz.  in  an 
oncoming  flow  Mach  number  M.  from  0  -  0.5.  The  drive  was  installed  in  the  FML 
indraft  wind  tunnel  at  .NASA  Ames  Research  Center.  Stroboscopic  schlieren  studies 
and  interferograms  as  well  as  holographic  interferometry  studies  were  conducted  for 
a  wide  range  of  flow  conditions,  amplitudes  and  frequencies.  The  results  show  that 
compressibility  effects  appear  at  .M  -  0.3.  that  a  dynamic  stall  vortex  forms  for 
all  Mach  numbers,  and  that  for  M  y  0.3.  the  dynamic  stall  angle  decreases  as  M 
increases.  On  the  otherhand.  increasing  the  degree  of  unsteadiness,  i.e.  increasing  the 
reduced  frequency,  delays  deep  stall  monotonically  at  all  Mach  numbers.  Amplitude 
has  a  dominant  effect,  and  dynamic  stall  occurs  at  lower  angles  of  attack  at  lower 
amplitudes,  yet  this  angle  of  attack  is  always  higher  than  the  static  stall  angle. 

Interferometry  studies  have  yielded  some  exciting  new  pictures  of  the  density 
field  around  the  leading  edge,  which  are  being  currently  analyzed. 
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2.  INTRODUCTION 

The  phenomenon  of  dynamic  stall  is  very  important  in  the  field  of  helicopter 
aerodynamics.  The  present  day  helicopter  flight  envelope  is  severely  restricted  be¬ 
cause  of  the  adverse  effects  that  are  produced  during  dynamic  stall,  notably  the 
excursions  in  the  pitching  moment  associated  with  the  streamwise  movement  of  the 
dynamic  stall  vortex,  with  the  result  that  its  beneficial  effects  remain  unharnessed. 
The  problem  of  dynamic  stall  is  well  known  and  relates  to  the  dynamic  delay  of  stall 
by  production  of  dynamic  lift,  when  an  airfoil  is  rapidly  pitched  past  the  static  stall 
angle.  Under  such  unsteady  motion  of  the  airfoil,  a  large  dynamic  stall  vortex  is  pro¬ 
duced  with  its  origin  near  its  leading  edge  and  converts  over  it.  The  presence  of  the 
vortex  contributes  significantly  to  the  dynamic  lift,  but  the  convection  of  the  vortex 
towards  the  trailing  edge  results  in  severe  dynamic  loads  on  both  the  pilot  and  the 
helicopter,  and  hence,  the  aircraft  controls  are  designed  to  restrict  the  vehicle  from 
living  into  those  condit  ions.  A  comprehensive  review  of  dynamic  stall  has  been  given 

(. 


airfoil  chord 

frequency  of  oscillation.  Hz 
reduced  frequency  = 
free  stream  Mach  number 
free  stream  velocity 
chordwise  distance 
angle  of  attack 
mean  angle  of  attack 
amplitude  of  oscillation 
circular  frequency,  radians/sec 


by  Carr'.  The  phenomenon  is  very  complicated  due  to  its  dependence  on  several  fac¬ 
tors  namely,  the  airfoil  shape,  the  leading  edge  geometry,  the  degree  of  unsteadiness 
(expressed  in  terms  of  the  reduced  frequency),  free  stream  Mach  number,  Reynolds 
number,  three  dimensionality,  the  state  of  the  boundary  layer,  the  amplitude  of  the 
airfoil  motion  and  mean  angle  of  attack,  etc.(  McCroskey'^).  The  process  of  dynamic 
stall  is  greatly  affected  ^v  compressibility.  Earlier  experiments'^  have  shown  that  the 
compressibility  may  even  ncj^a'^e  the  benefits  of  dynamic  stall.  However,  these  early 
studies  did  not  explore  the  problem  completely.  Some  computational  studies  (Fung 
and  Carr"',  Ekaterinaris'’ )  have  sho'vn  that  a  supersonic  flow  region  forms  over  the 
airfoil  surface  even  at  low  .Mach  numbers.  McCroskey  et  aP  have  shown  that  at  a  free 
stream  Mach  number  as  low  as  0.2.  the  flow  can  develop  high  suction  peaks  leading 
to  the  formation  of  a  local  supersonic  region.  The  consequent  appearence  of  shock 
waves  and  their  interaction  with  the  boundary  layer  could  have  a  dramatic  effect  on 
the  dynamic  stall  process.  It  is  very  important  to  understand  the  basic  physics  of  the 
process,  so  that  eventually  a  means  can  be  devised  for  its  control  and  exploitation  in 
real  systems.  This  report  describes  an  experimental  effort  that  was  carried  out  as  a 
first  step  in  this  direction. 


3.  FAC  ILITY.  TNSTRTTMENTATION  &  TECHNIQUE 

3.1.  Description  of  the  Drive  .System 

The  experiments  were  conducted  in  the  in-draft  wind  tunnel  of  the  Fluid  Me¬ 
chanics  Laboratory  (FML)  at  N.-\S.A  .\mes  Research  (Center  (.4RC).  It  is  one  of  the 
ongoing  dvnatnic  siall  research  projects  of  tlie  .Nav\-N.\SA  Joint  Institute  of  Aero- 
nautcis  between  the  Naval  Postgraduate  School  and  N.\SA  ARC. 

The  details  of  the  FML  iii-draft  wind  tunnel  are  given  in  Carr  and 
C'handrasekhara'  .  The  facility  is  one  of  a  complex  of  4  in-draft  wind  tunnels  con- 


nected  to  a  108  m'^  sec  (240,00  CFM).  9.000  hp  evacuation  compressor.  The  test 
section  size  is  25.4cmX  35cm.  X  100cm.  The  flow  in  the  tunnel  is  controlled  by  a 
variable  cross  section  dowmstream  diflPuser.  Its  throat  is  always  kept  choked  so  that 
no  disturbances  can  propogate  upstream  from  the  other  tunnels  or  the  compressor 
into  the  test  section. 

A  unique  mechansim  was  designed  and  built  to  produce  the  oscillatory  motion 
of  the  airfoil.  It  is  described  in  Ref.  7.  The  drive  system  is  loacted  on  top  of  the 
test  section.  The  test  section  windows  are  connected  to  the  drive  and  the  sinusoidal 
movement  of  the  windows  results  in  an  identical  movement  of  the  airfoil  which  is 
supported  by  the  windows.  The  airfoil  is  supported  by  pins  push  fitted  between  two 
2.54cm.  thick  optical  quality  glass  windows.  The  airfoil  supports  are  unique  in  that 
the  pins  carry  the  entire  load.  The  pins  are  smaller  than  the  local  airfoil  thickness 
and  hence  provide  complete  optical  access  to  the  airfoil  surface.  This  makes  detailed 
flow  studies  possible  even  at  the  surface. 

The  oscillating  drive  was  designed  to  meet  the  following  specifications: 

o  =  Oi,  -  a„,sjr?27r//  =  q,.  —  amSinuit 

0  <:  Q.j  <  15'’ 

2  '  1  <  10  ' 

0  ■_  /  <  100// r 
0  ::  <  0.5 

200.000  ■; /?c  -  10’ 

airfoil  chord  7.()2cni. 

The  airfoil  chord  and  the  flow  c  onditions  correspond  to  a  helicopter  in  forward  flight 
and  the  other  conditions  correspond  to  those  of  a  ith  scale  model  rotor,  whose  test 
results  are  directly  applicable  to  a  helicopter  rotor. 


Considerable  effort  was  expended  in  designing  a  system  capable  of  delivering 
the  specifications.  Special  care  was  needed  in  the  design  to  transfer  the  dynamic 
aerodynamic  loads  on  the  airfoil  to  the  glass  windows  that  support  it.  In  view  of  the 
material  properties  of  glass,  special  sleeves  of  delrin  with  wall  thickness  of  2.54/im 
had  to  be  custom  machined  and  inserted  between  the  glass  and  the  pins  that  support 
the  airfoil  to  cushion  the  direct  loads  on  the  glass.  The  drive  system  was  carefully 
designed  to  maintain  the  same  linkage  ratio  at  all  amplitudes  and  angles  of  attack 
so  that  the  variation  in  the  harmonic  motion  is  always  less  than  0.1%.  A  flywheel 
system  insures  the  constancy  of  the  harmonic  motion.  The  basic  load  pattern  on  the 
drive  is  cyclical  and  thus  fatigue  life  became  an  important  factor  to  be  considered. 
The  weakest  elements  in  the  system  are  the  airfoil  support  pins.  Other  similar  load 
bearing  components  include  the  bearing  pins  in  the  push  rods,  the  pins  connecting 
the  windows  to  the  push  rods.  The  drive  system  and  its  weakest  components  were 
all  designed  for  a  fatigue  life  of  10'  cycles  at  100  Hz  oscillation  frequency.  In  view  of 
the  complex  geometry  of  some  of  the  components  and  loads  on  them,  considerable 
effort  had  to  be  spent  to  develop  a  reliable  design.  The  system  is  driven  by  an  a.c. 
motor,  which  is  connected  to  the  drive  by  a  transmission  bell.  All  the  design  details 
re  documented  in  an  internal  stress  analysis  report  generated  for  the  NASA  ARC 
facility  approval  process.  The  facility  with  its  drive  system  is  referred  to  as  the 
Compressible  Dynamic  Stall  Facility,  (CDSF).  Fig.  1  shows  the  CDSF. 

3.2.  Phase  Locking  Instrumentation 

The  CDSF  is  instrumented  with  3  digital  encoders.  Of  these,  one  is  an  absolute 
position  encoder  providing  3600  counts  per  revolution  and  is  used  for  the  mean 
angle  of  attack  information.  The  other  two  are  incremental  position  encoders  with  a 
resolution  of  800  counts  per  revolution.  One  of  the  incremental  encoders  is  used  for 
obtaining  the  frequency  /  phase  angle  information.  The  other  could  be  used  for  the 
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instantaneous  angle  ol  attack.  But.  for  the  series  of  experiments  being  reported,  it 
was  not  used. 

A  triggering  circuit  was  specailly  built  for  obtaining  conditionally  sampled  data 
of  the  oscillating  airfoil  flow  field.  One  of  the  modules  in  this  circuit  is  a  compara¬ 
tor,  which  compares  the  desired  phase  angle  set  as  a  BCD  number  corresponding 
to  the  appropriate  encoder  count  with  the  changing  encoder  counts  streaming  into 
the  module.  When  a  match  occurs,  the  circuit  outputs  a  TTL  pulse,  which  in  turn 
triggers  a  light  source  for  the  schlieren  studies  or  a  laser  for  interferometry  studies. 
The  actual  phase  angle  at  which  this  occurs  is  immediately  frozen  on  the  display 
for  documentation.  .No  significant  delays  were  noticed  in  the  firing  of  the  schlieren 
light  source.  However,  when  the  laser  was  triggered  by  this  means,  the  finite  lasing 
time  resulted  in  the  actual  phase  angle  being  different  from  the  selected  number.  To 
properly  account  for  this,  the  circuit  was  modified  to  freeze  the  displays  only  when 
the  laser  flashes  -  by  using  a  very  sensitive  photo  diode  that  senses  the  laser  light 
flash.  In  view  of  the  different  oscillating  frequencies  at  which  the  airfoil  was  oscil¬ 
lated.  the  delay  constantly  changed.  Hence,  on-line  calibration  was  performed  for 
each  frequency  and  amplitude  of  oscillation  and  the  information  was  used  appropri¬ 
ately.  The  trigger  circuit  also  accomodates  a  display  that  indicates  the  frequency  of 
oscillation  and  outputs  this  information  for  input  to  a  computer. 

.Since  detailed  LD^  studies  are  also  planned,  additional  demands  were  placed 
on  the  phase  locking  circuitry  by  the  random  nature  of  the  LDV  data.  Since  the 
LD\  data  rate  is  dependent  on  particle  arrival  rate,  which  is  in  general  random  and 
depends  upon  the  local  flow,  a  means  to  read  the  instantaneous  phase  angles  arises 
each  time  an  LD\'  data  sample  is  validated.  This  requires  latching  circuits  to  freeze 
the  continuous  encoder  data  based  on  an  event  in  the  flow-  which  is  the  appearance 
of  co-incident  LDV  data  in  the  I  and  V  directions.  Fig.  2  presents  a  schematic  of 
the  method  followed.  The  LD^’  data  is  input  to  a  NASA  LDV  multiplexer,  to  which 
the  encoder  outputs  are  also  connected.  The  co-incidence  detection  pulse  (i.e.  the 
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data  ready  pulse)  from  the  multiplexer  is  used  to  freeze  the  encoder  data,  until  all 
the  data  is  completely  transferred  to  a  micro  VAX  II  computer  in  the  DMA  mode. 
As  soon  as  this  is  accomplished,  the  latches  are  released  for  fresh  data.  In  view 
of  the  general  paucity  of  the  LDV  data  and  the  time  it  takes  to  transfer  it  to  the 
computer,  very  high  speed  latches  are  needed  both  from  reasonable  experiment  time 
and  difficulty  considerations.  Such  components  have  been  used  in  this  circuit.  Fig. 
3  shows  a  partial  schematic  of  the  circuit  developed. 

3.3.  Details  of  Stroboscopic  Schlieren  Technique 

The  schlieren  system  used  is  standard  and  is  shown  in  Fig.  4.  It  is  also  described 
in  ref.  7  and  8.  The  phase  angle  at  which  the  schlieren  picture  was  desired  was  set 
as  the  corresponding  encoder  count  on  the  front  panel  switches  of  the  phase  locking 
and  triggering  circuit  as  a  BCD  number.  Pictures  were  then  taken  at  this  phase 
angle.  A  movie  of  the  flow'  was  also  obtained  for  the  case  of  M  =  0.3,  k  =  0.05, 
o  =  lO'"'  +  10‘''stnu;t.  To  do  this,  another  electronic  circuit  was  designed  to  control 
the  shutter  of  a  stop-action  camera  for  a  set  number  of  limes.  After  the  required 
number  of  frames  for  persistence  were  exposed,  a  different  phase  angle  was  selected 
and  the  process  repeated  through  a  complete  cycle  of  oscillation. 

3.4.  Stroboscopic  Dark  Center  Interferometry 

Dark  Center  Interferometry  is  a  technique  which  enables  real  time  interfero- 
grams  of  the  flow  to  be  obtained.  The  real  lime  aspect  of  the  technique  makes  it  a 
very  attractive  alternative  to  more  conventional  methods  such  as  holography.  The 
simpler  requirements  on  the  quality  of  the  optical  components  makes  it  a  preferred 
technique  over  the  traditional  Mach-Zehnder  interferometry.  However,  its  applica¬ 
tion  to  unsteady  flow  problems  presents  significant  challenges  both  in  terms  of  the 
requirements  on  the  light  source  and  on  the  recording  equipment. 


The  technique  uses  nearly  the  same  optical  configuration  as  a  schlieren  system 
and  is  shown  in  Fig.  5.  The  optical  components  on  the  transmitting  side  of  the 
system  provide  a  parallel  beam  of  light  which  has  been  expanded  to  fill  the  windows 
after  spatially  filtering  it.  The  beam  coming  out  of  the  test  section  is  focussed  to  a 
spot  on  a  transparent  coated  film  strip.  The  coating  on  the  film  when  exposed  to  the 
intense  energy  of  the  laser  beam  burns  to  create  a  dark  spot,  which  blocks  most  of 
the  light  through  it.  This  spot  is  created  with  no  flow  through  the  test  section.  When 
the  flow  is  turned  on.  the  deflection  of  the  beam  due  to  the  density  changes  causes 
it  to  be  focussed  around  the  spot  rather  than  on  the  spot  itself.  The  diffracted  light 
from  around  the  spot  and  the  little  light  passing  through  the  spot  itself  then  interfere 
resulting  in  a  system  of  fringes  in  real  time.  The  fringes  satisfy  all  the  relationships 
of  the  standard  interferometry  technique.  Thus,  a  distribution  of  the  flow  density 
field  results.  By  suitable  fringe  analysis  techniques,  the  flow  field  quantities  can  be 
recovered. 

When  used  for  the  oscillating  airfoil  dynamic  stall  flow  field  studies,  after  the 
creation  of  the  spot,  the  flow  was  established  and  with  the  airfoil  oscillating,  pictures 
were  obtained  in  a  manner  similar  to  the  stroboscopic  schlieren  technique.  As  already 
staled,  the  laser  was  triggered  at  the  desired  phase  angle  of  motion  taking  into 
account  the  appropriate  delays  and  the  displays  were  frozen  when  the  light  appeared. 


3.5.  Experimental  Conditions 


The  matrix  of  experimental  conditions  for  stroboscopic  flow  visualization  is 
shown  in  Table  1.  For  selected  conditions,  real-time  interferometry  and  LDV  studies 
were  also  carried  out.  The  results  are  discussed  in  the  next  section. 
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4.  RESULTS  AND  DISCUSSION 


The  information  presented  below  has  been  drawn  from  the  various  papers  that 
resulted  from  the  funding.  As  such,  they  are  only  briefly  discussed  and  the  original 
sources  properly  referenced.  In  cases  where  such  a  reference  is  still  not  available  in 
open  literature,  more  details  are  discussed. 

4.1.  Calibration 

The  calibration  of  the  facility  included  the  calibration  of  the  wind  tunnel  and  the 
encoders.  The  wind  tunnel  was  always  run  with  its  variable  cross  section  downstream 
throat  choked  to  prevent  any  downstream  disturbances  from  propagating  upstream 
into  the  test  section.  The  calibration  process  determined  the  range  of  the  throat 
heights  that  are  possible  when  a  stable  shock  forms  in  the  diverging  portion  of  the 
diffuser.  Further  details  can  be  found  in  Carr  and  Chandrasekhara^. 

Calibration  of  the  encoders  was  done  by  resting  the  airfoil  on  a  machined  block 
on  which  the  airfoil  surface  was  machined.  The  block  could  be  set  at  different  angles 
of  attack  by  using  precision  wedges.  The  procedure  consisted  of  recording  the  encoder 
counts  for  each  of  the  angles.  It  should  be  noted  that  each  angle  of  attack  occurs 
twice  during  a  cycle  and  the  encoder  was  calibrated  for  both  the  upstroke  and  the 
downstroke  of  the  airfoil  motion.  The  data  showed  that  the  variation  of  angle  of 
attack  followed  a  sine  wave.  This  can  also  be  seen  from  Fig.  6  wherein  the  actual 
angle  of  attack  variation  has  been  read  by  a  computer  and  plotted  for  a  frequency  of 
20  Hz.  and  for  the  motion  described  by  o  =  lO*^’  ^ 

4.2.  Stroboscopic  Schlieren  Flow’  Visualization 

Stroboscopic  schlieren  technique  reveals  the  instantaneous  density  gradient  pic¬ 
ture  of  the  flow,  without  any  history  effects  in  it.  Each  picture  is  a  snap-shot  of  the 


1.3 


flow  at  the  instant  the  picture  is  taken. 

Fig.  7  presents  a  sequence  of  schlieren  pictures  for  the  case  of  a  =  10®  +  10®sinw^ 
at  a  reduced  frequency  of  0.1  for  M  =  0.3.  It  should  be  noted  that  at  low  Mach 
numbers,  the  density  gradients  are  small  except  near  the  leading  edge.  However, 
since  the  system  used  was  very  sensitive  the  flow  details  could  be  photographed. 

The  oncoming  flow  reaches  stagnation  at  a  point  slightly  below  the  leading  edge, 
which  appears  as  a  dark  region  in  all  the  frames  in  the  figure  around  the  leading  edge. 
A  movie  taken  for  a  similar  condition  (at  k  =  0.05)  clearly  shows  the  downstream 
movement  of  the  stagnation  point  on  the  pressure  side  of  the  airfoil  as  it  pitches 
beyond  an  angle  of  attack  of  10®  during  the  upstroke.  Eventually  by  a  ss  13®  —  14®. 
it  settles  down  at  «3-5%  chord  location. 

As  the  airfoil  pitches  up.  the  flow  around  the  leading  edge  accelerates  and  the 
suction  peak  develops.  This  region  appears  as  a  bright  region,  which  terminates  in 
another  dark  region  on  the  upper  surface.  As  the  angle  of  attack  increases  further, 
the  dark  region  grows  and  starts  convecting  downstream.  Whereas  the  exact  origins 
of  the  dynamic  stall  vortex  could  not  be  clearly  determined,  it  is  strongly  believed  to 
be  in  this  region  downstream  of  the  suction  peak  where  all  the  vorticity  input  due  to 
the  unsteady  motion  of  the  surface  accumulates.  By  x/c  =  10-15%,  the  dynamic  stall 
vortex  becomes  clearly  visible  as  a  tightly  wound  structure,  for  example  at  o  =  15.2® 
in  Fig.  7.  At  q  =  15.9  ,  the  gradients  within  the  vortex  begin  to  clearly  appear  as 
a  sharp  transition  line  separates  the  bright  region  from  the  dark  region  of  positive 
density  gradient.  The  flow  is  still  bounded  by  the  leading  edge  shear  layer  enveloping 
the  vortex  and  the  downstream  flow  is  attached  as  can  be  seen  at  a  =  15.2®.  where 
the  boundary  layer  downstream  thickness  has  slightly  increased.  Further  increase 
in  the  angle  of  attack  causes  the  vortex  to  move  downstream,  while  simultaneously 
thickening  the  boundary  layer.  Eventually,  at  o  =  18.1®.  the  vortex  is  shed  and  deep 
stall  occurs.  The  airfoil  generates  lift  till  this  angle  of  attack,  which  is  substantially 
higher  than  the  static  stall  angle  of  12.4".  Beyond  the  deep  stall  angle,  the  flow  is 
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completely  separated  from  the  leading  edge,  three  dimensional,  and  highly  turbulent, 
with  a  large  wake.  As  the  airfoil  pitches  down,  reattachment  occurs  at  considerably 
below  the  static  stall  angle.  In  the  figure,  it  can  be  seen  that  even  at  q  =  8.4°.  the 
process  is  not  complete.  The  movie  also  shows  the  highly  turbulent  region  as  deep 
stall  occurs  and  as  the  airfoil  reaches  the  top  of  the  cycle. 

Similar  features  were  observed  at  most  Mach  numbers,  reduced  frequencies  and 
amplitudes.  Chandrasekhara  and  Carr**,  Carr  et  al°,  Chandrasekhara  and  Brydges^°, 
and  Brydges^^  provide  these  details. 

4.3.  Shock  Formation  over  the  Airfoil 

There  is  a  general  consensus  among  researchers  that  the  rapid  airfoil  movement 
generates  such  strong  suction  pressures  that  the  local  flow  becomes  supersonic  and 
eventually  the  flow  returns  to  the  subsonic  state  via  a  shock.  Computational  studies 
by  Fung  and  Carr**  have  pointed  to  shock  formation.  However,  it  has  remained  a 
postulate  till  now  and  has  never  been  directly  verified  experimentally.  Fig.  8a  shows 
the  first  such  documentation  of  a  shock  in  the  flow  at  M  =  0.45,  k  =  0.075  and 
o  =  10''  -1-  2^’sinujt  at  o  =  9.1".  The  shock  is  located  at  x/c  »  7%  and  is  about  1.5/X 
high.  Although  its  position  is  unsteady,  the  above  numbers  reasonably  quantify  the 
shock  features.  It  remains  on  the  surface  till  a  =  9.6".  when  it  breaks  down.  Fig.  8b 
clearly  shows  that  the  breakdown  occurs  as  a  A  shock,  indicating  the  laminar  nature 
of  the  flow.  Generally,  no  shock  induced  separation  could  be  found,  but  occassionally 
a  thin  shear  layer  seemed  to  separate  when  the  shock  moved  into  the  dark  region 
downstream.  At  the  condition  tested,  it  was  a  highly  intermittent  phenomenon  and 
could  not  be  properly  documented  due  to  the  synchronizing  problems.  However,  a 
visual  observation  showed  that  it  was  occurring  occassionally  and  it  is  schematically 
depicted  in  Fig.  8c.  It  is  hoped  that  further  experiments  with  a  high  speed  camera 
will  allow  capturing  of  this  event  properly. 

1.5 


4.4.  Effect  of  Mach  Number 


The  primary  emphasis  of  this  research  is  the  study  of  the  compressibility  effects 
on  dynamic  stall.  This  effect  was  assessed  in  two  ways.  Firstly,  the  dynamic  flow  field 
pictures  were  compared  at  different  Mach  numbers,  at  the  same  reduced  frequency, 
when  the  dynamic  stall  vortex  was  approximately  at  50%  chord  location.  Such 
a  qualitative  comparsion  is  shown  in  Fig.  9  (from  ref.  8).  It  can  be  seen  that 
for  M  <  0.3.  the  angle  of  attack  at  which  the  vortex  reaches  x/c  =  0.5  is  15.1*^' 
and  independent  of  Mach  number.  For  M  >  0.3,  the  corresponding  angle  of  attack 
decreases  progressively  from  13.8'  at  M  -  0.35  to  12.3*^  at  M  =  0.45,  indicating  strong 
effects  of  compressibility  beyond  M  =  0.3.  Thus.  M  =  0.3  can  be  considered  to  be  the 
delineating  Mach  number  for  compressibility  effects.  These  results  were  quantified 
by  following  the  dynamic  stall  vortex  center  location  over  the  airfoil  surface  with 
increasing  angle  of  attack.  The  results  are  shown  in  Fig.  10  for  a  reduced  frequency 
of  0.05.  The  airfoil  angle  of  attack  is  plotted  as  a  function  of  the  chordwise  location  of 
the  vortex  in  it.  It  is  clear  that  within  the  experimental  scatter  and  the  subjectivity 
of  such  an  approach,  the  vortex  behavior  appears  identical  till  M  =  0.3.  Beyond  this, 
dramatic  changes  can  be  seen  as  the  vortex  begins  to  appear  earlier  in  the  pitching 
cycle,  i.e.  at  lower  angles  of  attack  as  M  increases.  Thus,  deep  dynamic  stall  also 
occurs  at  lower  angles  of  attack  as  the  Mach  number  is  increased.  This  result  leads 
to  the  conclusion  that  compressibility  promotes  the  dynamic  stall  process. 

4.5.  Effect  of  Reduced  Frequency 

The  effect  of  increasing  the  reduced  frequency  (or  the  degree  of  unsteadiness) 
is  shown  in  Fig.  11  wherein  once  again  the  vortex  center  location  is  plotted  as  a 
function  of  angle  of  attack  for  different  Mach  numbers.  The  curves  for  larger  k  are 
higher  indicating  that  as  the  reduced  frequency  increases,  the  vortex  remains  closer 
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to  the  leading  edge  at  higher  angles  of  attack.  Thus,  the  deep  stall  angle  steadily 
increases  as  k  increases.  For  example,  at  M  =  0.3,  the  deep  stall  angle  is  15.9°,  and 
for  k  =  0.1,  it  is  18.1°.  This  result  is  true  even  at  higher  Mach  numbers,  and  in 
the  compressibility  range.  Even  though  it  was  seen  from  the  previous  section,  that 
increasing  Mach  number  promotes  stall,  increasing  k  has  the  opposite  effect  at  all 
Mach  numbers.  This  is  due  to  the  fact  that  the  vorticity  generated  by  the  airfoil 
motion  is  larger  at  higher  k  and  thus  is  available  for  the  vortex.  This  additional 
vorticity  enables  the  airfoil  to  produce  lift  till  higher  angles  of  attack  until  the  vortex 
is  shed.  Thus,  increasing  k  alleviates  dynamic  stall.  Ref.  8  and  9  provide  additional 
results. 

4.6.  Effect  of  Amplitude 

One  of  the  issues  still  largely  unknown  in  dynamic  stall  is  the  effect  of  the  am¬ 
plitude  of  unsteady  motion.  McCroskey°  states  that  the  effect  is  large.  To  assess  this 
effect  in  general  and  in  particular  when  compressibility  effects  are  present,  experi¬ 
ments  were  conducted  at  amplitudes  of  2°.  5"  and  10°.  Ref.  9  describes  the  results  of 
this  study  in  detail  and  the  salient  features  are  summarized  below.  Fig.  12  presents 
the  effects  at  different  Mach  number  and  reduced  frequencies.  In  all  these,  the  up¬ 
permost  curve  is  for  the  largest  amplitude  of  10°  and  the  lower  most  curve  is  for  the 
2°  amplitude.  It  is  evident  that  at  larger  amplitudes,  stall  is  delayed  to  higher  angles 
of  attack.  This  is  a  somewhat  surprising  result  because,  for  example,  at  =  5°,  for 
M  =  0.35.  the  dynamic  stall  angle  is  12. 9‘  ,  w’hereas  for  =  10°.  it  is  15.2°.  One 
expects  that  the  flow  at  the  lower  amplitude  would  be  still  attached  at  15.2"  because 
the  degree  of  unsteadiness  is  the  same  in  both  cases.  However,  on  further  analysis, 
it  becomes  clear  that  the  flow  is  subjected  to  a  much  larger  surface  acceleration  at 
the  larger  amplitudes.  Also,  the  local  flow  pressure  gradients  are  larger  for  the  larger 
amplitudes  because  of  the  motion  history,  for  motion  starting  from  the  same  mean 
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angle.  Thus,  the  total  amount  of  vorticity  generated  is  larger  at  the  larger  amplitudes 
leading  to  the  conclusion  that  such  flows  can  sustain  greater  pressure  gradients  and 
thus,  the  vortex  remains  over  the  surface.  The  result  appears  to  be  valid  both  in 
the  incompressible  and  compressible  flow  regimes.  It  is  also  interesting  to  point  out 
that  the  dynamic  stall  vortex  still  formed  for  all  cases  studied,  even  though  in  some 
cases,  the  maximum  angle  of  attack  reached  is  still  below  the  static  stall  angle.  Such 
a  case  is  shown  in  Fig.  12  where  the  vortex  has  just  formed  very  near  the  top  of  the 
cycle  for  the  case  of  q  =  10*^  -h  2^ sinujt  at  M  =  0.35. 

4.7.  Interferometry  Studies 

Fig.  13  presents  one  of  the  several  interferograms  obtained  for  the  oscillating 
airfoil.  The  experimental  conditions  for  this  are  M  =  0.45.  k  =  0.05,  a  =  10.32*^  and 
the  angle  of  attack  was  varied  as  a  =  10''  -i-  lO'^'stnujt.  The  stagnation  region  can  be 
seen  on  the  lower  surface  at  «  2  -3  %  chord  location.  In  the  accelerating  flow  from 
here,  the  presence  of  a  strong  suction  can  be  inferred  from  the  closely  spaced  fringes 
and  the  Mach  number  is  increasing  rapidly.  A  large  cell  of  accelerating  flow  forms  at 
x/c  ss  6  -  7%.  where  the  local  Mach  number  reaches  about  0.8  in  the  last  discernible 
fringe.  The  shear  layer  underneath  the  cell  shows  a  large  number  of  fringes  as  the 
outer  flow  merges  with  the  boundary  layer.  A  closer  examination  of  the  shear  layer 
reveals  the  curvature  of  the  fringes,  indicating  that  dynamic  stall  has  already  started. 

The  inset  in  the  figure  provides  additional  details  of  the  flow  near  the  leading 
edge.  The  large  number  of  fringes  in  this  region  once  again  indicates  that  the  density 
(and  pressure)  gradients  are  significant  here.  Currently,  software  and  algorithms  are 
being  developed  for  improving  the  fringe  quality  digitally  and  to  provide  automated 
analysis  so  that  these  details  can  be  recovered  from  the  large  number  of  interferograms 
obtained. 
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4.8.  LDV  Studies 


LDV  studies  have  begun  for  obtaining  the  velocity  data  of  the  streamwise  and 
vertical  components  of  velocity  in  the  dynamic  stall  flow  field.  In  addition  to  the 
complexities  described  in  the  section  on  instrumentation,  there  are  the  complications 
arising  from  the  requirements  of  proper  seeding  and  the  fact  that  close  to  the  airfoil, 
blockage  of  the  beams  occurs  during  part  of  the  oscillation  cycle  and  hence  no  data 
exists.  Further,  owing  to  the  fact  that  the  measurements  are  being  made  of  the 
unsteady  flow  field,  a  large  number  of  samples  have  to  be  collected  to  arrive  at  a 
meaningful  ensemble  average  at  each  of  the  desired  phase  angles.  Typically,  the 
data  shown  in  Fig.  14  were  obtained  from  10,000  samples  in  each  of  the  velocity 
channels  (along  with  the  phase  positions).  The  data  were  sorted  into  36  bins  by 
software  developed  specially  for  this  experiment  and  the  averages  computed  for  each 
bin.  The  software  package  is  very  comprehensive  and  incorporates  several  features 
for  data  validation,  binning,  curve  fitting  with  a  Fourier  series  of  variable  number 
of  terms  so  that  data  in  missing  bins  can  be  interpolated,  histogram  plotted  etc. 
Data  acquisition  is  ongoing  with  a  nominal  flow  grid  involving  20  midspan  vertical 
points  at  20  streamwise  stations  amounting  to  a  total  of  4  million  samples.  Given  the 
random  nature  of  LDV  data  and  the  fact  that  the  valid  data  (coincident  data  in  both 
the  U  and  V  channels)  rate  is  generally  poor,  ranging  from  a  few  to  a  few  hundred 
(400  at  best)  per  second,  the  experiments  are  very  time  consuming  and  difficult. 
Nevertheless,  the  effort  has  proceeded  reasonably  succesfully  as  can  be  seen  from 
Fig.  14  obtained  for  M  =  0.3,  k  =  0.05  and  o  =  10"  -  lO'  sinwf.  at  x/c  =  0.0,  y/c  = 
0.167,  z/c  =  0.0.  The  streamwise  veloctiy  varies  from  about  115  m/s  to  145  m/s  as 
the  phase  angle  from  0  degrees  (q  -  10")  to  ^  210  At  a  =  O",  i.e.  when  the  phase 
angle  is  90  degrees,  the  flow  velocity  reaches  the  flow  velocity  corresponding  to  M  = 
0.3  and  then  rises  to  its  maximum  value  as  the  airfoil  pitches  up.  The  velocity  drops 
rapidly  to  95m/ s  at  a  phase  angle  of  about  215  degrees,  which  corresponds  to  the 
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dynamic  stall  angle  of  about  15.9"  as  determined  from  the  schlieren  studies,  (see 
Table  2).  This  agreement  is  very  good.  The  velocity  recovers  to  the  0  degree  phase 
value  as  the  airfoil  returns  to  its  mean  angle  of  attack.  A  corresponding  variation  in 
the  vertical  velocity  can  also  be  seen  with  its  minimum  at  90*^  phase,  (q  —  0").  As  the 
stall  angle  approaches,  the  vertical  velocity  component  rises  dramatically  following 
the  separating  leading  edge  streamline.  Once  again,  as  the  airfoil  returns  to  its  mean 
angle  of  attack,  the  flow  reattaches  and  the  vertical  veloity  drops  to  the  zero  phase 
angle  value.  Measurements  at  other  locations  in  the  flow  field  show  that  the  flow 
is  affected  at  large  distances  from  the  airfoil  indicating  the  strong  interactions  that 
exist  between  the  inner  viscous  and  the  outer  inviscid  regions.  These  are  the  first  set 
of  dynamic  stall  velocity  data  obtained  with  non-inirusive  optical  instrumentation 
and  it  is  hoped  that  in  a  very  short  while,  the  substantial  velocity  data  base  that  is 
currently  being  generated  will  lead  to  some  new  ideas  about  the  flow  field. 

5.  CONCLUSIONS 

The  study  has  led  to  the  following  major  conclusions: 

1.  A  dynamic  stall  vortex  is  present  at  all  .Mach  numbers. 

2.  Compressibility  has  a  significant  effect  on  the  dynamic  stall  process.  The 
delineating  Mach  number  between  the  incompressible  and  compressible  flow  regimes 
is  0.3.  Stall  is  promoted  by  increasing  the  Mach  number  above  this  value. 

3.  A  shock  is  present  on  the  airfoil  in  the  region  0  <  ^  <  1%  for  certain  flow 
conditions. 

4.  Increasing  reduced  frequency  alleviates  the  stall  process. 

5.  Increasing  the  amplitude  results  in  the  dynamic  stall  being  delayed  to  higher 
angles  of  attack.  At  all  amplitudes.  Mach  number  and  reduced  frequency  have  similar 
effects. 

6.  Interferometry  studies  indicate  strong  gradients  within  the  dynamic  stall 
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vortex.  The  local  Mach  number  has  been  estimated  to  reach  a  value  as  high  as  0.8 
at  a  free  stream  Mach  number  of  0.45. 


7.  The  leading  edge  region  supports  a  number  of  fringes  in  the  interferograms, 
pointing  to  the  large  suction  peaks  and  accelerations  there. 

8.  LDV  studies  offering  the  first  velocity  data  in  the  dynamic  stall  flow  field 
under  compressibility  conditions  have  been  obtained. 

Further  analysis  of  all  the  results  is  under  way  presently. 
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TABLE  1:  EXPERIMENTAL  CONDITIONS 
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Steady  flow 
10  deg.  ampl itude 
5  deg.  amplitude 
2  deg.  amplitude 


TABLE  2:  VORTEX  RELEASE  ANGLES  OF  ATTACK 
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First  row:  10  deg.  amplitude 
Second  row:  5  deg.  amplitude 
Third  row:  2  deg.  amplitude 
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SIDE  PANEL  WINDOW 


Fig.  1.  Schematic  of  the  Compressible  Dynamic  Stall  Facility 
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.  Partial  Schematic  of  the  Phase  Locking  and  Triggering  Circuit 
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Fig.  4.  Schematic  of  the  CDSF  Schlieren  and  LDV  Instrumentation 
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ARRANGEMENT  ABOVE 


Fig.  5.  Schematic  of  the  Optical  Arrangement  of  the  Dark  Center  Interferometer 
System  for  the  CDSF 
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Fig.  6.  Representative  Calibration  Data  for  the  Oscillating  Airfoil  Drive:  q 
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Fig.  7.  Stroboscopic  Schlieren  Photographs  of  the  Compressibility  Effects  on  Dy¬ 
namic  Stall  of  an  Oscillating  Airfoil:  M  =  0.3,  k  —  0.1.  a  =  10"  -  10^5tn(^t) 
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Fig.  9.  Effect  of  Compressibility  on  Dynamic  Stall  of  an  Oscillating  Airfoil:  k  =  0.05 


11.  Kffccts  of  Kodured  Frequency  on  Dynamic  Stall:  (a)  M  =  0.25.  (b)  M 
r  i  M  ().:i5.  (d)  M  0.45.  o  10'  - 


Fig.  i:;.  ai-'l'inu'  Intcrfcrograrn  of  the  Dynamic  Stall  Flow  Over  an  Oscillating 
Airfoil:  M  k  O.Q').  n  lO''  -  \0"sinu.'t  =  10.32”,  (Inset  shows  details 


around  the  leading  edge) 


